Polarized States and Domain Walls in Spinor Bose-Einstein Condensates 
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We study spin-polarized states and their stability in ant i- ferromagnetic states of spinor (F = 1) 
quasi-one-dimensional Bose-Einstein condensates. Using analytical approximations and numerical 
methods, we find various types of polarized states, including: patterns of the Thomas-Fermi type; 
structures with a pulse-shape in one component inducing a hole in the other components; states 
with holes in all three components; and domain walls. A Bogoliubov-de Gennes analysis reveals 
that families of these states contain intervals of a weak oscillatory instability, except for the domain 
walls, which are always stable. The development of the instabilities is examined by means of direct 
numerical simulations. 



I. INTRODUCTION 

The development of the far-off-resonant optical tech- 
niques for trapping ultracold atomic gases has opened 
new directions in the studies of Bose-Einstein conden- 
sates (BECs). Atoms can be confined regardless of 
their hyperfine (spin) state, thus avoiding freezing the 
atom's spin degree of freedom, making it available for 
the study of interesting spin dynamics [1]. One of the 
major achievements in this direction was the experimen- 
tal creation of spinor BECs ^2|, j3|] . The spinor conden- 
sate formed by atoms with spin F is described by a 
(2F + l)-component macroscopic (mean- field) wave func- 
tion, which gives rise to various phenomena that are 
not present in single-component BECs, including forma- 
tion of spin domains [4], spin textures and multi- 
component (vectorial) solitons of bright [y, S S], dark 
0, gap and brigfit-dark [ll| types. 

Generally, the dynamics of the spinor F = 1 BEC is 
spin-mixing [12]. However, there also exist non-mixing, 
or spin-polarized states of the system, which are sta- 
tionary solutions of the corresponding system of Gross- 
Pitaevskii (GP) equations [13, 14]. The stability of such 
polarized states is different for the two distinct types of 
the F = 1 condensates, namely the ferromagnetic (FM), 
such as in ^'^Rb, and polar (alias anti-ferromagnetic^ 
AFM), such as ^^Na, ones, where the spin-dependent in- 
teractions are, respectively, attractive and repulsive. Ac- 
cordingly, as demonstrated in Ref. [15] (see also Ref. [7]), 
spin-polarized states are modulationally stable/unstable 
in the AFM/FM condensates. 

In this work, we focus on AFM spinor condensates, and 
study, in particular, spin-polarized states of the spinor 
BEC of ^^Na atoms. Assuming that this spinor sys- 
tem is confined in a strongly anisotropic trap, we first 
present the respective system of three coupled quasi one- 
dimensional (ID) GP equations. Then, employing the 
so-called single- mode approximation 0, [l5| , we use an- 
alytical and numerical methods to find the spin-polarized 



states of the system and study their stability via the 
Bogoliubov-de Gennes (BdG) equations (i.e., the lin- 
earization of the GP equations). The simplest possi- 
ble form of these states are the stable Thomas-Fermi 
(TF) configurations (for each hyperfine component). We 
also present other spin-polarized states, including ones 
in which one component maintains a pulse-like shape, 
inducing a hole in the other components, and structures 
with holes in all the three components. All these states 
feature instability regions too, but with values of the 
normalized instability growth rate up to ~ 10 De- 
velopment of the oscillatory instabilities is examined by 
means of direct simulations. It is found that the struc- 
tures with a hole in one component get weakly deformed, 
while the states with holes in all the three components 
suffer stronger deformations. If the three components are 
initially spatially separated, being confined in different 
harmonic traps, but are put in a single trap afterwards, 
spin domain-wall (DW) patterns are formed. A family of 
the DW solutions exists and is fully stable if, for a fixed 
value of the trap's strength, the chemical potential (or 
the number of atoms) exceeds a certain critical value. 

The paper is organized as follows: Section II presents 
the model. Sections III is dealing with TF states. In Sec- 
tions IV and V we examine structures with one and multi- 
ple holes, respectively (including their stability). Finally, 
conclusions are presented in Section VI. 



II. THE MODEL AND SETUP 

At sufficiently low temperatures, and in the framework 
of the mean-field approach, a spinor BEC with F = 1 
is described by a vectorial order parameter, ^(r, t) = 
[^_i(r,t),^o(r,t),^+i(r,t)]^, with the different fields 
corresponding to three values of the vertical component 
of the spin, rriF = — 1,0,+1. Assuming that this conden- 
sate is loaded into a strongly anisotropic trap, with hold- 
ing frequencies Ux <C cj^, we assume, as usual, that the 
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wave functions are separable, ^o,±i = i^o,±i{x)i^^{y, z), 
where the transverse components ip±{y,z) represent the 
ground state of the respective harmonic oscihator. Then, 
fohowing the standard approach [16] of averaging the cou- 
pled 3D GP equations for the three components in the 
transverse plane z), leads to the system of coupled ID 
equations for the longitudinal components of the wave 
functions (see Refs. [6, 7, 9, 1^,0]): 



ZW±1 = i^si^/'±l+4'°^(|V'±l|' + IV'0|'-|V'Tl| 



.(ID) .,2 



„(1D), 



24'^V-i^5V'+i, 



+ iV'+ir)^o 



(1) 

(2) 



where star denotes complex conjugate, and Hs] 



.(ID) 



ntot is the spin- 



independent part of the Hamiltonian, with ntot = 
l^/^-ip -h |'0oP + IV^+iP being the total density (m is 
the atomic mass). The nonlinearity coefficients have 
an effectively ID form, namely Cq^^^ = co/27ra^ and 

c>2^^ = C2/27ra^, where a± = y^h/muj± is the trans- 
verse harmonic-oscillator length which determines the 
size of the transverse ground state. Coupling constants 
Co and C2 account for, respectively, the mean- field spin- 
independent and spin-dependent binary interactions be- 
tween identical spin-1 bosons, 

{co, C2} = {4:7Th^/3m) {{a2 + ao), (a2 - ao)} , (3) 

where ao and a2 are the s-wave scattering lengths in com- 
bined symmetric collision channels of total spin / = 0, 2. 
The spinor F = 1 condensate with C2 < and C2 > is, 
respectively, of the FM _and AFM types, as, respectively, 
in ^^Rh and ^^Na [il,ll2|]. 

Measuring time, length and density in units of 



h/cQ^^riQ^ h/ymcQ^^riQ and no (here, no is the peak 
density), we cast Eqs. ([I])-(j2j) in the dimensionless form. 



±1 



+ ^ror^i. (4) 

idt^o = i/siV^O + ^(|V^-l|' + |V^+l|')V^O 

+ 25^^iro^^u (5) 

where i^si = -(l/2)9^ + (l/2)(^^^x^+ntot, the normalized 
trap's strength is 



2(ao + 2a2)no \u;± 
and the parameter 5 is given by 



I2 ~ O'O 



(1^) ao + 2a2 



(6) 



(7) 



{S < and ^ > correspond, respectively, to the FM 
and AFM spinor condensate). For spin-1 ^^Rb and ^^Na 



atoms, this parameter is (5 = —4.66 x 10~^ [3 ^ = 
3.14 X 10~^ [19], respectively. Thus, in these physically 
relevant cases, (5 is a fixed small parameter of the system. 

Spin-polarized states of the system, characterized by 
a constant population of each spin component, can be 
found upon substitution 



i^j = ynj{x)eyip{iOj - ijiijt) , j = -l,0,+l, (8) 

where nj and Oj are densities and phases of the compo- 
nents, and /ij are their chemical potentials. Substituting 
this in Eqs. 0|) and ([5|), it is readily found that conditions 
for the existence of spin-polarized states are 



2/io 
AO 



2^0- 



= or TT. 



(9) 
(10) 



Below, we assume that the chemical potentials of the 
components are equal, ja = jii-i = jhq = /i+i. For anal- 
ysis of the stability of a stationary spin-polarized state, 
^sps(^) = [V^-i(x), Vio(x), Vi+i(x)]^, we win perform lin- 
earization around the unperturbed spin-polarized state, 
assuming a perturbed solution. 
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^pj{x, t) = i}j{x) + e (uj{x)e~^^* + v*{x)e^'^*^ 

(11) 

where uj and vj represent infinitesimal perturbations 
with eigenvalues A = + ^A^. Then, the solution of 
the ensuing linear-stability problem (i.e., the BdG equa- 
tions) for A and associated eigenfunctions Uj , v* provides 
complete information about the stability of the under- 
lying stationary state, ^sps- Whenever it is unstable, 
we will also examine its evolution through direct simula- 
tions of GP equations (|4]) and ([5]), using a finite difference 
scheme in space and a fourth-order Runge-Kutta time in- 
tegrator. Typically, in the simulations the unstable state 
is initially perturbed by a uniformly distributed random 
perturbation of relative amplitude 5 x 10~^. 

To estimate relevant physical parameters, we assume 
the spinor condensate of ^^Na atoms with peak ID den- 
sity no = 10^ m~^, confined in the harmonic trap with 
frequencies lj± = 27r x 230 Hz and cj^^ = 27r x 13 Hz. 
In this case, the normalized trap strength is = 0.1 
(this value is kept fixed throughout this work), while the 
number of atoms, depends on chemical potential fi 
and the particular form of spin-polarized states. Typi- 
cally, we take /i in interval 1 < /i < 5, which corresponds 
to values of N in the range of 3.5 x 10^ < N < 3.5 x 10"^ 
atoms. For instance, at /i = 2 the number of atoms is 
N ~ 10^; in this case, the normalized time and space 
units correspond, respectively, to 1.2 ms and 1.83 /im. 



III. THOMAS-FERMI SPIN-POLARIZED 
STATES 

The simplest spin-polarized states can be found in the 
framework of the single-mode approximation jTH, ITEI ]. 
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In anticipatio n, sug gested by the approximation that 
components Uj {x) may be close to eigenmodes of a 
single effective potential, induced by the combination 
of the trap and nonlinearity, we introduce the ansatz 
nj{x) = qjn{x)^ where coefficients qj are the populations 
of each spin component in the steady state, related by 
the normalization condition, Q-i -\- qq -\- g+i = 1. Then, 
Eqs. dH)-© lead to the following system: 



L^S 1 



— {sqo - 2^q-iq^i) ] n 
Q±i 



V^ = 0, (12) 



L^S{l-qo^2s^q-iq+i)n = 0, (13) 



where L = -ii- {l/2)dl + {l/2)^'^x'^ + n and 5 = ±1 for 
= or = TT, respectively. It is readily observed 
that the three equations ([T2|), ([13]) are tantamount to a 
single one, viz. 



pn ] ^Jn = 0, 



(14) 



supplemented by relations 

p=l + (5, qo = 2^qZTq^ for A^ = 0, (15) 

p = 1, q-i — q-^i for AO = tt. (16) 

In uniform space {Q = 0), Eq. (p!4|) leads to the con- 
stant density, n = fi/p. As shown in Refs. [3, El, [lH, 
such constant solutions to Eqs. (jl])-© are modulation- 
ally stable only in the AFM phase, with S > 0. Be- 
low, we will only consider the case of the AFM conden- 
sate (in particular, bosonic sodium atoms with spin 1 
and S = 3.14 X 10~^). In the presence of a sufficiently 
weak trap, Eq. (p!4|) can be solved approximately employ- 
ing the TF approximation [20]: neglecting the kinetic- 
energy term (~ d^^/n) in Eq. (p!4|) . we find density pro- 
files of the three spin components: in the region where 
/i > (1/2)1]V, 



(17) 



and rij = elsewhere. Obviously, all three components of 
the TF solution have the same spatial width, or Thomas- 
Fermi radius, Rtf = \f2\ilVt. 

In numerical simulations, we used a fixed-point al- 
gorithm (Newton-Raphson method) to find exact spin- 
polarized solutions to Eqs. (|4|)-([5]), with profiles close to 
those produced by the TF approximation, as given by 
Eq. (pT|) . In particular, we used, as an initial guess, three 
identical profiles of the form 



?/^j(x) = \J n{x) exp(z6>j). 



(18) 



with n{x) = /i — (l/2)r^^x^ and A^ = or tt. Then, keep- 
ing the trap's strength, fixed, we varied the chemical 
potential /i, and the numerical solution converged to sta- 
ble spin-polarized states, which were indeed close to the 
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FIG. 1: (Color online) Examples of two stable spin-polarized 
states for = (left panel) and A^ = tt (right panel), 
both obtained for Q — 0.1 and fi — 2. The wavefunctions 
ip±i are identical and are depicted by the solid line, while the 
wavefunction ipo is depicted by the dashed line. In the left 
panel, q-i = = 0.5 and = 1, while in the right panel 
q-i — q-^i — 0.25 and qo — 0.5. 



approximate one given by Eq. (pT|) . Two typical exam- 
ples are shown in the top panels of Fig. [T] for both cases, 
A^ = (left panel) and A^ = tt (right panel), with 
= 0.1 and ji = 2. The numerically determined states 
are very close to their TF-predicted counterparts, with 
q-i = g+i = 0.5 and qo = 1 {AO = 0, left panel) and 
q-i = q-^i = 0.25 and qo = 0.5 (A^ = tt, right panel). 



IV. SPIN-POLARIZED STATES WITH HOLES 

Apart from the TF spin-polarized states, there exist 
other ones, which feature holes in some of the compo- 
nents. An example of such states is shown in Fig. [21 As 
seen in the two top panels of this figure (both a stable and 
an unstable state are shown — see below), the ipo compo- 
nent is concentrated in the form of a pulse located at the 
trap's center, while the components feature a large 
hole at the same spot. This arrangement is explained 
by the fact that the interaction between components is 
repulsive, hence a peak (hole) in i/jq (V^±i) induces a hole 
(peak) in tlj±i (V^o)- The norm N of each component 
is shown, as a function of chemical potential /i, in the 
third-row panel of Fig. [2l The set of the linear stabil- 
ity eigenvalues for this state is shown in the second-row 
panels of Fig. [21 for two different values of fi: The state 
with /i = 3 (top-right) is stable, as all the eigenvalues are 
imaginary, while the state with ja = 2 (top-left panel) is 
unstable. In fact, all such unstable states are destabi- 
lized by a Hamiltonian Hopf bifurcation, which leads to 
a quartet of eigenvalues with nonzero real parts. The in- 
stability interval is 1.81 < ja < 2.15, with the maximum 
instability growth rate Max(Ar) ~ 1.3 x 10~^, found at 
/i 2 (see bottom panel of Fig. [2]). 

A similar state with the roles of ip±i and ipo exchanged, 
i.e., the i/jq component featuring the hole, and ones 
concentrated in narrow pulses, were also found. More- 
over, such states were found with either tjjo or having 
the opposite sign (i.e., for A^ = tt). The results are not 
shown here, as the (in)stability of these states is qualita- 
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FIG. 2: (Color online) The spin-polarized state with a hole in 
each of the ip±i components, and a pulse- like ipo component, 
for Q = 0.1 and = 0. The two top left and right panels 
show unstable and stable states, with ii — 2 and /j = 3, 
respectively; solid and dashed lines depict the components 
V^ii and ipo. The two panels in the second row display the 
spectral planes (Ar, A^) of the (in) stability eigenvalues for the 
same states. Note that the instability (at /i = 2) is of the 
oscillatory type, accounted for by a quartet of eigenvalues 
with nonzero real parts. The panel in the third row shows the 
normalized number of atoms (norm), N, of each component 
as a function of chemical potential /x; solid and dashed lines 
show the norms of components V^ii and -00, respectively. The 
bottom panel shows Ar as a function of /i, which reveals the 
instability window at 1.81 < fi < 2.15, with a maximum 



instability growth rate (Xr)^ 



1.3 X lO"-" at /i 2. The 



latter value corresponds to the unstable state shown in the 
top-left panel. 



tively the same as in the above case. 

The evolution of unstable states is exemplified in 
Fig. [3l which displays results of direct simulations of 
Eqs. (|4|) and ([5]) for the unstable state with = 
and /i = 2, that was presented in Fig. [2] (top- left panel). 
In Fig. [3l contour plots of densities of the components 
of the solution are displayed as a function of time (the 
densities of the and i/j-i components are identical). 




FIG. 3: (Color online) Top panels: Contour plots of densities 
of a solution shown in Fig. [2] subject to the oscillatory instabil- 
ity. The left and right panels display, respectively, identical 
densities of 'ip±i and that of ipo. Because of the extremely 
small growth rate of the instability, it manifests itself only 
for t > 4000. Bottom panels: The respective wavefunction 
profiles at t = (solid lines) and t = 10000 (dashed lines); as 
above, the left and right panels show, respectively, the wave- 
functions ip±i and '00- 



It is clearly observed that the predicted oscillatory in- 
stability sets in at a very large time {t ^ 4000, which 
corresponds to t ~ 5 seconds in physical units); this is a 
consequence of the extremely small growth rate of the in- 
stability. Eventually, the system settles down to a steady 
state, which has many similarities with the initial one. 
In particular, as seen in the top panels of Fig. [3l after 
t ^ 8000 the pulse in the i/jq component broadens and 
its amplitude is accordingly decreased, while the hole in 
the ip±i components becomes shallower. It is also noted 
that the three components develop a similar structure in 
their tails, as shown in the bottom panels of Fig. [3l 

Apart from the states considered above, it is also pos- 
sible to find spin-polarized states which feature, e.g., one 
hole in each of the components, and two holes in 
?/^o- Examples of such a state are shown in the top pan- 
els of Fig. [4] (the left one, for /i = 2, is stable, while 
the right one, for /i = 3 is unstable — see below). As 
seen in this figure, one may consider ip-i and as 
built of two overlapping pulses, which induce two holes 
in the ?/^o component due to the repulsive interatomic 
interactions. Results of the stability analysis for these 
states are shown in Fig. [H In this case, two quartets 
of eigenvalues with nonzero real parts are found in the 
spectral plane (see the second row right panel in Fig. [4|) 
in the interval 2.58 < /i < 3.22 (the bottom panel in 
Fig. [4]). The respective largest instability growth rate 

is (Ar)jnax ^ ^ ^ of same 

order of magnitude as in the previous case. The devel- 
opment of the instability was studied, as above, in direct 
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FIG. 4: (Color online) Same as Fig.[2]but for a state with one 
hole in each of the 2p±i components and two holes in the ipo 
component. In this case, the instability is manifested by two 
eigenvalue quartets the largest of which has a nonzero real 
part in the interval 2.58 < fi < 3.22; the maximum instability 
growth rate is (K)^^^ ^ 1.8x 10~^ for /x ^ 2.9. The unstable 
state shown in the top-right panel corresponds to /i = 3. 



simulations, starting with initial conditions in the form 
of a perturbed solution pertaining to /i = 3. The result 
is shown in Fig. O in terms of the evolution of identi- 
cal densities of the iIj±i components, and the density of 
tpQ. Again, the instability manifests itself at large times 
{t > 3500, which corresponds to t > 4.2 seconds in phys- 
ical units), but in this case the final result is a strong 
oscillatory deformation of the three spin components (af- 
ter t ~ 5500), contrary to what is the case in Fig. [3] (the 
establishment of a new spin-separated state). 

Similar states with one hole in tpo and two holes in 
each of the iIj±i components, as well their counterparts 
with = TT, have also been found. They are not shown 
here; as in the previous case, the (in) stability of these 
additional states is similar to that reported in Fig. [H 
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FIG. 5: (Color online) Same as Fig. O for the unstable state 
shown in the top-right panel of Fig. 2] (pertaining to /i = 3) . 
The instability manifests itself at large times (t > 3500) and 
results in strong oscillatory deformation of the clouds. 
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FIG. 6: (Color online) Initialization of the system when three 
different traps of the same strength but of different location 
are used: the TF clouds for V^o and are confined in 
harmonic traps of strength Q = 0.1 located, respectively, at 
X = — 10, X = (dashed lines) and x = 10 (solid line). The 
traps at X = —10 and x = are turned off and this spinor 
configuration, confined solely in the trap located at x = 10, 
is introduced in the fixed-point algorithm. 



V. DOMAIN WALLS 

In the above sections, we reported the spin-polarized 
states in which all three spin components were spatially 
overlapping, since they were confined to the same poten- 
tial trap. However, it is also possible to use three different 
traps, each confining a different component, to initially 
separate them, and then allow the system to evolve in the 
presence of one of these traps (i.e., turning off the other 
two). In this section, we present spin-polarized states, 
including DW structures, obtained in this way. 

First, we describe the initialization of the system. We 
assume that the three TF-shaped components are ini- 
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FIG. 8: (Color online) Evolution of the domain-wall structure 
shown in the top-left panel of Fig. [G] to which a random per- 
turbation was added. Shown in the left and right panels are 
spatiotemporal contour plots of the densities of components 
ip±i (identical ones) and ipo, respectively. 
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FIG. 7: (Color online) Top-left panel: the wave functions of 
the components (V^+i and ip-i are identical) in a stationary 
state found from the initial configuration prepared as shown 
in Fig. [6l The resulting spin-polarized state has the form of 
a domain- wall like structure between the ipo and iIj±i com- 
ponents. The parameters are Q = 0.1 and /x = 2. Top-right 
panel shows the wavefunctions of the domain- wall state found 
Sit fi — 1.43. Bottom panel: the norm of each component in 
the above-mentioned domain-wall structure versus the chemi- 
cal potential (the dependences for = or tt are identical). 
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FIG. 9: (Color online) Same as Fig.[8]but for the state shown 
in the top-right panel of Fig. [6] The left, middle and right 
panels show, respectively, the densities of ip-i, ipo and ip+i. 
Noteworthy is a stationary dark- solit on- like structure, located 
at X ~ 13 in the component. 



tially loaded into three different traps, Vj{x)^ of the same 
strength, centered at different positions: 

Vj{x) = ^n\x-jAxf, j = -l,0,+l. (19) 

We choose Ax = (i.e.. Ax = 10 for Q = 0.1), which 
implies initially overlapping TF configurations, see Fig.[6l 
After constructing this state, we turn off the traps V-i{x) 
and Vo{x)^ retaining only the rightmost one, y^i(x), 
which now acts on all the three components, and feed 
the initial configuration state into the fixed-point algo- 
rithm, to find new spin-polarized states. Other possibil- 
ities, such as turning off potentials V±i and keeping Vq, 
arranging the three components in a different way, etc., 
eventually lead to retrieving the spin-polarized states pre- 
sented in the previous sections, while the approach out- 
lined above [keeping V+i(x) and switching V-i{x) and 
Vo(^) off] generates new DW patterns, which are dis- 
played in Fig. [7J and could not be obtained otherwise. 
(The asymmetry of the procedure is instrumental in gen- 
erating the new states). 

The most interesting spin-polarized DW states found 
following this procedure correspond to values of the 
chemical potential /i > 1.43 (or norm N > 5400 for 
Q = 0.1); for smaller values of /i we typically find struc- 
tures of the Thomas-Fermi type. Two examples, one for 
/i > 1.43, and another corresponding to /i = 1.43, are 
shown in Fig. [3 In the former one, the tpo component 



(which has the larger norm) is centered to the right of 
the midpoint of the remaining trap (x+i = 10), while 
the identical iIj±i components are pushed to the left, due 
to the repulsion from i/jq, with a DW created between 
ip±i and t/jQ. Figure [3 displays an example of such a DW 
structure for /i = 2 and Q = 0.1. 

The state found at the above-mentioned value, /i = 
1.43, is shown in the top-right panel of Fig. [71 for AO = 0. 
In this state, the shape of the tpo component is similar 
to that considered above, while the iIj±i components are 
not identical, in contrast to the above example. In this 
case the i/jq and iIj±i components overlap over a wider 
spatial region. Note that i/j-i changes its sign at x 13, 
featuring a structure resembling the waveform of a dark 
soliton embedded in a bright one |2l|. 

The stability of the DW states was also investigated 
by means of the BdG equations. It was concluded that 
there are no eigenvalues with a real part in the interval 
1.43 < /i < 5, or, equivalently, 5400 < N < 35000 for 
Q = 0.1 (not shown here in detail). Thus, the DWs 
are stable in this region. Verification of the stability, 
performed by direct simulations of Eqs. (|3|) and ([5]), is 
illustrated in Figs. [8] and [9l for /i = 2 and /i = 1.43, 
respectively. It is obvious that these states are indeed 
stable, for very large times exceeding t = 10000 (i.e., 12 
seconds in physical units). 
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VI. CONCLUSIONS 

In this work, we have studied spin-polarized states in 
anti- ferromagnetic spinor {F = 1) Bose-Einstein conden- 
sates. In particular, our analysis applies to a quasi- ID 
spinor condensate of sodium atoms. The considerations 
were based on analytical calculations and numerical com- 
putations of the coupled Gross-Pitaevskii equations for 
this setting. 

Assuming that all three hyperfine (spin) components 
are confined in the same harmonic trap, we have found 
various types of spin-polarized states and examined their 
stability. The first family consists of Thomas- Fermi con- 
figurations, considered analytically in the framework of 
the single-mode approximation (which assumes the sim- 
ilarity of the spatial profiles of all the components). 
Within their existence region, these states were found to 
be stable. Also identified were more complex patterns, 
which include states composed of pulse-like structures in 
one component, which induce holes in the other com- 
ponents, and states with holes in all three components. 
These states feature windows of weak instability. The de- 
velopment of the instability was investigated by means of 
direct numerical simulations, which demonstrate that it 
manifests itself at very long times, and results in a weak 
deformation of the states with a single hole in some of 
the components, and a stronger one in the states with 
holes in all components. 

Fully stable families of spin-polarized states develop 



from initial configurations which are initially separated 
(by means of three different traps) components. These 
states form domain-wall structures between the compo- 
nents, at values of the chemical potential above a critical 
value. At the critical value (corresponding to a certain 
norm, for a fixed trap's strength), we have found another 
spin-polarized state (including a dark-soliton element), 
in which all the components partly overlap. 

It would clearly be interesting to investigate the ex- 
istence and stability of higher-dimensional counterparts 
of the ID spin-polarized states found in this work. An- 
other relevant question for further analysis is whether 
such spinor condensates support stable topological ob- 
jects, such as dark solitons or vortices. Work in these 
directions is in progress. 
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